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We investigate turbulence in dilute polymer solutions when polymers are strongly stretched by 
the flow. We establish power-law spectrum of velocity, which is not associated with a flux of a 
conserved quantity, in two cases. The first case is the elastic waves range of high Reynolds number 
turbulence of polymer solutions above the coil-stretch transition. The second case is the elastic 
turbulence, where chaotic flow is excited due to elastic instabilities at small Reynolds numbers. 
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INTRODUCTION 



In this paper we continue theoretical investigation of 
turbulence in dilute polymer solutions, started in [Q,^. 
As opposed to Newtonian fluids, such solutions possess 
additional macroscopic degrees of freedom related to the 
elasticity of the polymer molecules. Relaxation times of 
elastic stresses can be comparable with time-scales of the 
flow which means that the relation between the stress 
and the velocity gradient is non-local in time and, con- 
sequently, in space. It is a striking property of dilute 
polymer solutions that minute amounts of polymer can 
significantly modify the flow. Probably, the most famous 
example is the drag reduction phenomenon. Addition of 
long-chain polymers in concentrations as small as, say, 
10~^ by weight, can induce a substantial reduction of 
the drag force needed to push a turbulent fluid through 
a pipe ^1^. Another example is the elastic turbulence 
1^-^, which is a chaotic flow, excited in the dilute poly- 
mer solutions at low Reynolds numbers. 

The reason why small amounts of polymer can sig- 
nificantly modify properties of the fluid is flexibility of 
polymer molecules. At equilibrium a polymer molecule 
coils up into a spongy ball of a radius Rq. The value of 
Rq depends on the number of monomers in the molecule, 
which is usually very large. For a dilute solution with the 
concentration of the molecules, n, satisfying uRq ^ 1, an 
influence of equilibrium size molecules on the hydrody- 
namic properties of the fluid can be neglected. When 
placed in an inhomogeneous flow, such a molecule is de- 
formed into an elongated structure, which can be char- 
acterized by its end-to-end distance R. If the number 
of monomers in a typical polymer molecule is large, the 
elongation R can be much larger than Rq. The influence 
of the molecules on the flow increases with their elonga- 
tion and may become substantial when R^ Rq. 

Deformation of a polymer molecule is determined by 
two processes, stretching by the velocity gradients and re- 
laxation due to elasticity of the molecule. To understand 
how the molecule resists the deformation by the flow, one 



should consider its relaxation. Recent experiments with 
DNA molecules indicate that the relaxation is linear in 
the wide region of scales Rq <^ R <^ -Rmax, where i?max 
is the maximum molecule extension |^. In the case of 
polymers, theoretical arguments and numerics presented 
in jl^ support the linear relaxation. These results can 
be understood if we assume that at i? ^ i?o the role of 
excluded volume and hydrodynamic interactions between 
the monomers are negligible. Then the random walk ar- 
guments suggest that the entropy of polymer molecules 
is quadratic in R in the range i?o ^ -R ^ Rmax imply- 
ing linear relaxation. Whether the polymers are excited 
by the flow is determined by the softest relaxation mode 
that corresponds to the dynamics of the elongation R. 
In the absence of stretching, the relaxation of R is de- 
scribed by the equation dtR = —R/t, where r is a re- 
laxation time, which is expected to be i?-independent at 
Rq <^ R Rmax- If the end-to-end distance R is of the 
order of the maximum extension, t starts to depend on 
R and the dynamics of the molecule becomes nonlinear 
Iprl . Possible statistical consequences of the non-linearity 
have been investigated in 

The behavior of the molecule in an inhomogeneous 
steady flow depends on the value of the Weissenberg 
number, Wi, defined as the product of the characteris- 
tic velocity gradient and r. When a polymer molecule 
is placed in a flow, smooth at the scale R, the veloc- 
ity difference between the end-points is proportional to 
R multiplied by the characteristic value of the velocity 
gradient. At Wi ^ 1 the relaxation is fast as compared 
to the stretching time and the polymer always relaxes to 
the equilibrium size, Rq. The behavior of the polymer 
at Wi }^ 1 depends on the geometry of the flow. For 
purely elongational flows the molecule gets aligned along 
the principal stretching direction. If the velocity gradient 
is larger than the inverse relaxation time, i.e. Wi ^ 1, 
the elastic response becomes too slow in comparison with 
the stretching and the molecule gets substantially elon- 
gated The sharp transition from the coiled state 
to the strongly extended state is called the coil-stretch 
transition. Rotation can suppress the transition and even 
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damp it completely since the molecule does not always 
point in the stretching direction (see, e.g., For ex- 

ample, no coil-stretch transition occurs in the case of the 
shear flow, which is a particular combination of the clon- 
gational and rotational flows. 

In contrast to the steady polymer molecule, 

moving in a random flow, alternately enters regions of 
high and low stretching. As the intensity of the flow 
increases, the effect of the stretching becomes more pro- 
nounced. One can generally assert existence of the coil- 
stretch transition in this case. This has been first demon- 
strated by Lumley jlj], who considered the situation 
where the characteristic time of variations of the velocity 
gradient is much larger than the inverse of the character- 
istic value of the gradient. He showed that if the ampli- 
tude of the velocity gradient fluctuations is large enough, 
the expectation value of grows with time, which sig- 
nifles the coil-stretch transition. We have shown in [nU2] 
that the coil-stretch transition occurs in any random flow 
and established a general criterion for the transition. In 
particular, the transition occurs in the situation where 
the time of velocity gradient variation is of the order of 
the inverse of its characteristic value, which is likely to 
be the case for real flows. The coil-stretch transition in 
random flows is controlled by the parameter Air, where 
Ai is the average logarithmic divergence rate of nearby 
Lagrangian trajectories, to be referred to as the principal 
Lyapunov exponent (which is positive for an incompress- 
ible flow ||l^,^). The molecules are weakly stretched 
if AiT < 1 and strongly stretched otherwise. Therefore 
for random flows the parameter Air plays the role of the 
Weissenberg number. 

As it is well known (see, e.g., p7| , p^ ), turbulent flows 
in Newtonian fluids consist of chaotic eddies from a wide 
interval of scales, rj < r < L, where L is the integral scale 
(where the flow is excited) and ij is the viscous scale. The 
energy pumped at the scale L cascades down to the scale 
ry, where it dissipates. The size of the polymer molecules 
is usually much smaller than the viscous scale. Viscosity 
makes the flow smooth at scales r < rj, i.e. the velocity 
can be approximated by linear profiles at these scales. 
Therefore, if i? < 77, then the stretching of molecules 
is determined by the velocity gradient, which is random 
in the turbulent fiow. The Lyapunov exponent can be 
estimated as the characteristic value of the velocity gra- 
dient, which is determined by the eddies at the viscous 
scale rj. As the Reynolds number grows, the velocity gra- 
dient increases, and so does Air. At some value Rbc of 
the Reynolds number the product Air reaches the value 
1 and the coil-stretch transition occurs. 

Several mechanisms can limit the polymer stretching 
above the coil-stretch transition. The first one is the 
internal non-linearity of the elasticity of the polymer 
molecules. If this mechanism dominates, then above the 
transition the molecules are stretched up to the elonga- 
tion of the order of i?max- An alternative mechanism has 
been proposed by Tabor and de Gennes [|l9|. It is based 
on the assumption that the elongation of the polymer 



molecules R becomes larger than the viscous length of 
turbulence, ry, where the elastic force always wins over 
the stretching at a certain value of the elongation. Be- 
low, we assume that i? ^ 77, which seems to be rea- 
sonable for typical polymer solutions. Another limiting 
mechanism is the back reaction of the polymers on the 
flow. It is caused by the collective contribution of coher- 
ently deformed polymer molecules into the stress tensor. 
This elastic part of the stress grows with the molecule 
elongation. When it becomes of the order of the vis- 
cous stresses existing in the flow, the polymers modify 
the flow around them suppressing the stretching. As a 
result, a dynamic equilibrium is realized at a character- 
istic elongation, i?back- The total polymer stress is pro- 
portional to Tii?^, so that i?back depends on the polymer 
concentration n. We assume that the concentration is 
large enough for the value of -Rback to be much smaller 
than i?max. Probably, the condition i?back <C -Rmax is 
necessary for existence of a stationary state, because the 
polymer molecules, stretched up to i?max, are intensively 
destroyed by the flow. 

Above the coil-stretch transition the back reaction 
modifies the small-scale properties of turbulent flows, 
which leads to the emergence of a new scale, > rj, 
where energy dissipates mainly due to polymer relax- 
ation. The scale plays the role of a new dissipation 
scale. Large-scale eddies with the sizes r > do not 
excite elastic degrees of freedom so the usual inertial en- 
ergy cascade is realized at these scales. At Re ^ Rcc 
there appears a new region of scales, 77* < r < r*, where 
elastic waves can propagate [Q, which are analogous to 
the Alfven waves in magnetic hydrodynamics. At the 
scale 77* viscosity becomes essential, leading to the strong 
damping of the elastic waves. 

In this work we investigate the velocity spectrum in 
the elastic wave range fcr* ^ 1 (where k is the wave vec- 
tor). We show that the spectrum obeys a power law. The 
ideas of the analysis go back to the works of Townsend 
po| and Batchelor |2^]. They recognized that fluctua- 
tions with scales, smaller than the smoothness scale of 
the flow (77 in the case of usual turbulence and in our 
case), evolve in the linear velocity proflles. For the pas- 
sive scalar at large Prandtl numbers Batchelor derived 
the spectrum with the power-law | pl| , originating 
from the exponential character of stretching in the linear 
flow. Formally, it is explained by zero scaling dimension 
of the advection term that implies scale invariance. This 
property is not broken by a linear decay term where the 
power-law spectrum still holds, as it was shown for a lin- 
early decaying passive scalar by Chertkov ||2^. Here we 
introduce a consistent theoretical scheme for the descrip- 
tion of the small-scale fluctuations and show that though 
the dynamics of the small-scale fluctuations is more com- 
plicated, than in the case of the passive scalar (advection 
and linear decay accompanied by stretching and waves), 
the power-law in the spectrum still holds. The wave oscil- 
lations break the scale-invariance, but their influence on 
the energy balance is reduced to forcing the equipartition 
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of kinetic and elastic energies of small-scale fluctuations. 
The power law spectrum terminates at A: ~ (z/t)~^/^, 
where viscosity overcomes stretching. The power-law in- 
terval widens as Re grows. Let us stress that this law is 
not related to a flux of any conserved quantity. 

Another situation, where a power-law spectrum of the 
small-scale fluctuations is observed, is the elastic turbu- 
lence, realized in low- Reynolds polymer solutions, if the 
Weissenberg number Wi is large enough It was 

shown experimentally in |^-^ , that the coil-stretch tran- 
sition, occurring at increasing Wi, leads to a chaotic flow 
even though Re is small. Its existence is due to hydrody- 
namic instabilities caused by the presence of the elas- 
tic stresses. The velocity spectrum is observed to be 
power-like in a wide range of scales in this case 
We demonstrate that, in contrast to the usual hydro- 
dynamic turbulence, in the elastic turbulence the power 
velocity spectrum is not related to the energy cascade, 
since the main energy dissipation occurs at the largest 
scales. The mechanism, leading to this power spectrum 
is, again, similar to the linearly decaying passive scalar 
problem. 

The structure of the paper is as follows. In Section | 
we introduce a system of equations describing the coupled 
dynamics of inertial and elastic degrees of freedom. This 
system is similar to the system of equations describing 
the magnetohydrodynamics (MHD) [p3| with the impor- 
tant difference of a linear decay term in the equation on 
the "magnetic held" . In Section |l| we present results, 
concerning statistics of a passive scalar with a constant 
damping, embedded in a random flow. It is a prototype 
for the subsequent consideration. Section [II is devoted 



to the description of the principal properties of the large 
Re turbulence in the presence of polymers and to the 
derivation of the power-law spectrum in the elastic dissi- 
pation range. In Section IV we establish the power-law 
spectrum for the elastic turbulence. In Conclusion we 
summarize our results and discuss possible implications 
of our work for other subjects. Appendix is devoted to 
some details of the Lagrangian statistics. 



I. BASIC RELATIONS 

We study dynamics of dilute polymer solutions at 
scales much larger than the inter-molecular distance 
where the polymer solution can be regarded as a contin- 
uous medium and described by macroscopic fields. Char- 
acteristic times of considered processes are regarded to be 
comparable with the polymer relaxation time r. In this 
case, besides the usual hydrodynamic degrees of freedom, 
one has to take into account degrees of freedom, related 
to the polymer elasticity. These degrees of freedom can 
be described in terms of the polymer stress tensor Q . 

We assume that the flow can be treated as incompress- 
ible, that is V • D = 0, where v is the velocity of the flow. 
This is justified provided processes at a given scale are 



slow in comparison with sound oscillations at the same 
scale. Then the velocity dynamics can be described in 
terms of the following equation 0| 



(1.1) 



which is a generalization of the Navier-Stokes equation 
to the case of viscoelastic fluids. Here P is the pressure, 
J/ is the kinematic viscosity of the solvent, g is the fluid 
mass density, and Ily is the polymer contribution into 
the stress tensor per unit mass. 

The equation (|l . l| ) has to be supplemented by an equa- 
tion for the polymer stress tensor 11^ ||l^. We assume 
the following equation 



(1.2) 



where r is the polymer relaxation time, and Hq is related 
to the thermal fluctuations of the polymer conformations 
. Let us briefly repeat applicability conditions of the 
equation (|l.2| ) , discussed in Q . Linearity of the decay 
term in ( |1.2D assumes R <^ i?max, where R is the typ- 
ical polymer molecule size. The equation (1.2) implies 
that there is a single mode related to the polymer defor- 
mations, which is an idealization. A polymer molecule 
has a lot of deformational degrees of freedom, that have 
different relaxation times. A number of such degrees of 
freedom was observed in experiments with DNA ||9| . Nev- 
ertheless, in the turbulent flows, only the mode with the 
largest relaxation time is strongly excited, whereas other 
modes are excited at most weakly. Thus, Eq. (1.2) should 
be treated as the equation related to the principal mode. 

The conce ntratio n of the polymer molecules n enters 
the system (1.1 1.2) only via Hq, IIq oc n and implicitly 
via the assumption R <C Rmax since 11 cx nR^. If n is in- 



homogeneous, then the system of equations (1.1,1.2) has 
to be supplemented by the equation for the concentration 
dtn + vVn — (we neglect small diffusivity of polymer 
molecules) . In this paper we consider the case when poly- 
mer molecules are strongly ext ended. Then 11 ^ IIq, and 
the term with IIq in Eq. (L2) can be discarded. In this 
case any explicit dependence on the concentration of the 
polymer molecules n drops from the system of equations 
(1.1 L^ ). Therefore the dynamics of the polymer solu- 
tions with different values of n is identical in this regime 
as long as n is large enough for the condition R <C -Rmax 
to be satisfied. 



A. Lagrangian Description 

One can establish some properties of the polymer stress 
tensor 11, using the Lagrangian description of a fiuid. It 
is based on the notion of fluid particles trajectories (La- 
grangian trajectories) x{t,r), which are determined by 
the relations 



dtx = v{t,x), x{to,r)=r. 



(1.3) 
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The poi nt r plays the role of a Lagrangian marker. If IIo 
in Eq. (L2) is neglected, then it is possible to write its 
solution as 

n(<, x) = W{t, to)n(io, r)W^{t, to)e-2(*-*«)/" , (1.4) 

where the superscript T denotes a transposed matrix. 
Here W is the Lagrangian mapping matrix determined 
by the relations 



dtW{tM) = ^{t)W{t,to). 



(1.5) 
(1.6) 



Above a is the tensor of the velocity derivatives along the 
Lagrangian trajectory x(t) that includes the strain ten- 
sor and local rotations. The incompressibility condition 
V ■ w = is formulat ed i n terms of cr as tr a — 0. Then a 
consequence of Eq. (L5) is det W = 1. 

The matrix W describes deformations of infinitesimal 
fluid volumes. For example, the separation, 5x, between 
two close fluid particles, moving along the Lagrangian 
trajectory x{t), evolves according to 



5x{t) = W{t,t')5x{t') . 



(1.7) 



Therefore Wij{t,tQ,r) — dxjjt r)/dri. Now it is easy to 



understand a meaning of Eq. ( L4 ) . The polymers are ad- 
vected along the Lagrangian trajectories being stretched 
by the velocity gradient and relaxing to their equilibrium 
shape due to the polymer elasticity. 

We now briefly describe statistical properties of the 
matrix W, details can be found in Appendix. We rep- 
resent the matrix as W = MAN, where M and N are 
orthogonal matrices, while A is a diagonal matrix. At 
times much larger than the velocity gradients correlation 
time T(j the main eigenvalue exp(pi) of A becomes much 
larger than the rest, under the condition that the set of 
the Lyapunov exponents is non-degenerate. If a statis- 
tically steady state is realized, then the observa tion time 
is arbitrarily large, and we conclude from Eq. (1.4) that 
the matrix 11 has to be uniaxial 



Ilifc — BiBi. 



(1.8) 



as it was noted in g. This conclusion is almost self- 
evident once one goes back to the derivation of Eq. (1.2), 
recognizing that at i? 3> -Ro the contribution of ther- 
mal fluctuations into 11 is negligible, so that 11^^ cx RiRj 
holds. We observe that the vector B characterizes the 
direction and the strength of the coherent molecule elon- 
gations weighted by their contribution into the stress ten- 
sor. Note, that B is deflned up to sign, in analogy with 
the director in nematic liquid crystals. It follows from 
(O) and (O) that 



B{t, x) = exph(i - to)/T]W{t, to)B{to, r) 



(1.9) 



If t — to ^ T then W in this relation can be estimated as 



There are some modifications of the W statistics with 
respect to a Newtonian fluid, that are i mpo sed by the 
above relations. As it follows from Eq. (1.9), stationar- 
ity of the B statistics implies that pi{t) — t/r has a sta- 
tionary distribution. In particular, we conclude that the 
principal Lyapunov exponent Ai = limt^oc Pi/t of the 
flow is equal to 1/r, independently of the Reynolds num- 
ber. This means that above the coil-stretch transition 
the characteristic value of the velocity gradient is fixed 
at the scale 1/r. The above behavior is contrasted to the 
Newtonian fluids for which Ai grows with increasing Re 
and fluctuations of pi — t/r grow with time. The absence 
of the growth of the fluctuations is related to anticorrela- 
tions in the temporal dynamics of the component dii{t) 
of ffij along B. These anticorrelations show themselves 
in the equality / dt((a-ii(0)(Tii(i))) = (double brackets 
designate an irreducible correlation function) and origi- 
nate in the special interaction of the incrtial and elastic 
degrees of freedom, explained in more detail in Section 

m. 



B. Dynamic Equations for Dilute Polymer Solutions 



It is convenient to rewrite the equations (1.1,1.2) in 
terms of the vector B thus getting rid of extra degrees of 
f reed om. Substituting the decomposition (1.8) into Eq. 
(1.2), one obtains 



dtB, + {vV)B, = BkVkV^ - B./t . 



(1.10) 



This equation is similar to the one satisfied by the mag- 
netic field in MHD j2^, with the constant damping in- 
stead of the magnetic resistivity. The resemblance is 
made even stronger by notin g that B has to be solenoidal. 
Indeed, it follows from Eq. ( 1.1C| ) that 

dtV ■B + {vV)V -B^-V -B/t. (1.11) 

Therefore V • B monotonically decays becoming zero in 
the (statistically) ste ady state. Substituting the expres- 
sion ( pT^ ) into Eq. (1.1), and taking into account the 
constraint V ■ B ~ one obtains 

dtv + {vV)v = vV'^v - g-^VP + {BV)B . (1.12) 



Now the analogy of the system (1.10,1.12) with the sys- 
tem describing MHD at zero magnetic resistivity is 
almost comp lete. The only difference is in the damping 
term in Eq. ( |l.lC| ). 

The energy density per unit mass is given by the sum of 
kinetic /2 and elastic B^/2 t erms . The energy balance 
equation, following from Eqs. ( |1.10| ,1.12), is 

[dt+vV) (t.V2 + SV2 



{B\I){B-v) 



-vVi{vViv) 



(1.13) 



The energy dissipation is due to the viscous a nd th e poly- 
mer relaxation terms. Other terms in Eq. (1.13) repre- 
sent energy fluxes (in real space), they can be written as 
full divergences due to the constraints V • = = V • JB. 
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II. MECHANISM OF SCALE-INVARIANCE: 
PASSIVE SCALAR WITH LINEAR DAMPING 

Before investigating statistical properties o f the poly - 
mer solutions, described by the equations ( 1.1C , 1.12 ), 
we present statistical properties of a passive scalar with 
constant damping at scales smaller than the smoothness 
scale of the flow, examined in This simple case en- 
ables one to recognize the origin of a power spectrum for 
passively advected fields. 

The equation for the passive scalar 9 in the considered 
case is 



dte + v\je = -e/T 



(2.1) 



where r is the passive scalar decay time and 4>is a. forcing 
term needed to maintain the stationary state. It is as- 
sumed to be concentrated at a finite range of wavevectors 
near kf. We have omitted the diffusive term which can 
be neglected in comparison with the constant damping 
(in some region of scales) provided the diffusion coeffi- 
cient is small enough. Note that the constant damping 
of the passive scalar leads to a well-defined steady statis- 
tics even in the absence of the passive scalar diffusion. 

A flow is smooth on scales smaller, than the velocity 
gradients correlation length /. The smoothness means 
that a velocity difference between two points can be ap- 
proximated by a linear profile 



5vi 



TijSrj , 



(2.2) 



where 6rj is the separation between the points and a 
is a function of time. Obviously, (7^ — VjU^. For the 
usual turbulent velocity the correlation length I is equal 
to the viscous scale, I = 77. The smoothness of the turbu- 
lent velocity at scales less than 77 was first exploited by 
Batchelor [£1|, who considered statistical properties of a 
passive scalar at these scales. 

The linearity of the velocity difference leads to a power 
law for the passive scalar spectrum E(k), which is de- 
fined as {9{k)9{k')) = {2TT)^S{k + k')E{k), where 9{k) is 
Fourier transform of 9(r) with the wave vector fc, angular 
brackets designate averaging over the statistics of v, and 
we assume homogeneity and isotropy of the statistics. In- 
deed, the time of the energy transfer from k to 2k at kl ^ 
1 is scale-independent due to the linearity of the velocity 
profile. On the other hand, during the spectral trans- 
fer time the energy decay is also fc-independent (since 
the damping term is scale- independent). As a result the 
spectral fimction E(k) satisfies a relation E{2k) — CE{k) 
(with a constant C < 1). The solution of this equation 
is a power-law E{k) cx /c^" with 2"C — 1. 

Now we put the above consideration into a more rig- 
orous frame. We consider the passive scalar spectrum 
E{k) at fcZ ^ 1. The evolution of wave packets with 
such wavevectors is determined by the velocity gradient 
a. Let us consider the evolution during a time io a-nd 
express 9{to) via 0(0). A value of 9{to) near a point ri is 



determined by an evolution of 9 in the vicinity of the La- 
grangian trajectory a;(i,ri). To examine this evolution, 
one may perfo rm the Taylor expansion of the velocity 
V in Eq. (2.1) up to the first order in r — a; since the 
homogeneous advection does not affect equal-time cor- 
relation functions due to the Galilean invariance. Then 
one obtains 



dt9+[u + (T-{r~x)]\I9 = 



(2.3) 



Here u{t) = v{t, x) and a — a{t, x) are the velocity and 
the velocity gradients matrix along the Lagrangian tra- 
jectory X. Fourier transform 9k of the field 9 measured 
in the moving frame 9{t, r) = 9{t, r + x) satisfies 



d,9k^[ka--^ 



9k = - + <f>k exp [ik ■ x] . (2.4) 

T 



Further we confine ourselves to wavevectors k ^ kf, that 
is much larger than those on which the pumping (jj is sup- 
ported. In this case 9^ is determined by the convection 
from smaller wavevectors and the forcing term can be ne- 
glected. Then the equation (|]^) can be solved explicitly, 
and we find 



9{t,k) 



-t/r 



9(0, kW) 



(2.5) 



where W = W (t, ) is the Lagrangian mapping matrix, 
see Subsection [A. Returning to the real space, we ob- 
tain 



9{t,r + x) = e 



dfc 

(2^ 



.^k[r+W^m0(^Q^J^^r) ^ (2.6) 



The above formula is valid for |M^"^r| ^ I, the condition 
means that the passive scalar coming to a point r + ri 
was all the time in the l-vicinity of x allowing the Taylor 
expansion for the velocity. 

Let us now consider the pair correlation function of the 
passive scalar f{r) = {9{to,ri)9{to,ri + r)), defined as 
the spatial average over ri. We assume that the average 
(9) is zero (which can always be achieved by a shift of 
9 by a co nstant). The product of the fields is given by 
Eq. ( |2.6| ) [remind that a;(to) — fi] and depends on ri 
via the argument of W. The average over space (over 
ri) is equivalent to the average over the velocity statis- 
tics, or over the velocity gradients statistics along the La- 
grangian trajectories. If Ai^o S> 1 then the average over 
the interval < t < to and negative times can be done 
independently. Indeed, a{t) has a Lagrangian correlation 
time Aj~^. Thus velocity at negative times is correlated 
with a only at |t| ^ X^^ while W{to) is not sensitive to 
the value of a there, due to Aj^^ <C to (see Appendix). 
Therefore we can write 



fir) 



-2to/T 



dfc 

(2^ 



eM^kr){E[kW{to)]), (2.7) 



where E{k) is the spectrum function introduced above. 
Noting that E{k) equals the Fourier transform of the 
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pair-correlation function we obtain the following station- 
arity condition for the spectrum (we substitute to by t) 



E{k) = (exp(-2t/r)£;(A;VF)) 



(2., 



where W ~ VF(i, 0). The equation ( ^.^ ) is applicab le a t 
fcZ ^ 1, as follo ws f rom the condition r ^ Hn Eq. (2.6). 

The relation ( |2.8| ) has a simple meaning. The wavevec- 
tors of small-scale fluctuations of the passive scalar evolve 
according to k{t) = k{Q)W-^{t) as was shown by Kraich- 
nan Thus the energy of a fluctuation with the 

wavevector k is equal to its energy time t ago at the 
wavevector kW{t, 0) multiplied by the decaying factor 
exp(— 2i/T). Note that we could equall y we ll start di- 
rectly from ( p.4| ) to derive the equation ( p.S] ). One can 
show that in the spatially homogeneous situation one can 
always introduce such equation for the investigation of 
the spectrum at kl ^ 1. 

The equation (2^) is the quantification of the heuris- 
tic arguments given in the beginning of the section taking 
into account that the energy transfer time is by itself a 
random quantity. Its solution is a powe r law E{k) cx k~°' . 
Substituting the expression into Eq. ( |2.8| ), one gets the 
relation 



exp(2t/T) = {\kW/k\-°') , 



(2.9) 



which determines the exponent a. At Xit 3> 1 the mo- 
ments of |fcW^/fc| behave exponentially with time. Indeed, 
they are roughly equal to the product of Ait indepen- 
dent identically distributed random variables. Besides at 
these times the moments are independent of k/k due to 
the isotropization of W{to) described in Appendix. As a 
result the above equation has a uniqu e ph ysical solution 
examined in more detail in Appendix A3, where the in- 



equality a > 3 is established. The inequality has simple 
meaning that the spectrum has to decay faster than the 
Batchelor spectrum k~^ (fc~^ in the spherical normaliza- 
tion) holding at infinite r. 

Let us now extend the above results. The power law 
spectrum persists, even if the relaxation time r is k- 
dependent, but scales as zero power of k, that is if r 
depends on the direction of k only. This dependence 
can be regular (which makes the spectrum anisotropic) 
or random. Another remark is that addition of an oscil- 
lating term (with cjfe) into the equation for 0k, 



ka 



d_ 
dk 



9k/T-iuok{t)ek 



(2.10) 



does not change its spectrum even though the oscillating 
term breaks the scale-invariance. Indeed, let us pass from 
9 to 0, which is 9 — exp{iipk)9, with the phase, satisfying 



dtfk 



— tOk ■ 



(2.11) 



Then for 9 we return to the equation (2.6), that leads to 
the power spectrum. It remains to note that the spec- 
tra of 9 and 9 coincide. In other words, oscillating terms 



conserve energy and are, consequently, irrelevant for the 
energy balance. 

Below we generalize the simple picture, presented in 
this section, to the polymer solutions. 



III. HIGH REYNOLDS FLOWS 

Here we consider turbulence in dilute polymer solu- 
tions, when the Reynolds number exceeds the critical 
value RCc, corresponding to the coil-stretch transition. 
Then the polymer molecules are strongly elongated. Two 
different cases are possible, depending on the concentra- 
tion of the polymer molecules n. If it is very small, the 
elastic stresses are small in comparison with the viscous 
stresses. Then the polymers are stretched to their max- 
imal elongation, i?max, and the properties of the fluid 
do not differ significantly from those of the pure sol- 
vent. Below we consider the second, more interesting, 
case, when the concentration of the polymers n is large 
enough, so that elastic stresses can be larger than the 
viscous stresses. Then the polymer back reaction sub- 
stantially modifies the fiow. 

Whereas in the pure solvent typical velocity gradients 
grow unlimited as the Reynolds number increases, in 
polymer solutions above the coil-stretch transition the 
balance of inertial and elastic degrees of freedom fixes the 
characteristic value of the velocity gradient at 1/r. In- 
deed, if the instantaneous velocity gradient exceeds l/r, 
it extends the polymers, so that the elastic stress grows 
and damps the gradient. On the other hand, if the ve- 
locity gradient is smaller than l/r, the molecules con- 
tract and their influence on the flow diminishes. Then 
the velocity gradients tend to grow up to the value char- 
acteristic of the pure solvent, which is larger than l/r. 
Thus the velocity gradients fluctuate near l/r, that ex- 
plains the statistically steady state realized above the 
coil-stretch transition. Let us derive the condition for 
the existence of this steady state, related to the exis- 
tence of the maximal size -Rmax of the polymer molecules. 
In the vicinity of the coil-stretch transition \7v ~ l/r 
so th at i?back ^ vVv ~ v/t as it follows from Eq. 
(1.12). This leads to the condition for the existence of 



the back reaction regime -B^^x ^ ^/'''j where iJ^^x 
the maximal value of the elastic stress tensor achieved at 
R ^ i?max- Using estimates for the microscopic parame- 
ters, proposed in Ref. ||2^, one can rewrite this condition 
as n ^ (i?o-Rmax)~^- Below we will assume that there 
exists an interval in Re such that Bmax exceeds the value 
of B prescribed by the flow. The latter increases as Re 
grows so that the condition will break down at certain 
Re. After this happens either polymer degradation oc- 
curs or polymers start to behave as rigid bodies with size 
-Rmax- In the latter case the fluid becomes Newtonian 
again with renormalized viscosity. 

We assume Vt <^ L, where V is the velocity at the 
turbulence integral scale. Then the gradient related to 
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the large eddies is smaller than t~^. Therefore, the large 
eddies do not excite polymers, which means that the elas- 
tic stress tensor is not correlated at these scales. Since 
only coherent excitations of the elastic stress tensor can 
influence the flow, we conclude that the elasticity is neg- 
ligible for large eddies. The interaction of inertial and 
elastic degrees of freedom becomes essential at a scale 
r*, where velocity gradients are of the order of l/r. Here 
the energy starts to dissipate due to the polymer relax- 
ation, that is the inertial cascade terminates at r ~ r*. 
Since the velocity gradients fluctuate near the value l/r, 
reached at r* , at r < r, velocity difference scales linearly 
with the distance that is r^, is the smoothness scale of the 
flow. Near the coil-stretch transition characteristic veloc- 
ity gradient is determined by the viscous scale and is of 
the order of l/r, hence r* ~ r]. As the Reynolds number 
increases, velocity fluctuations increase, so that the scale 
r, grows which is very different from the Newtonian flu- 
ids where the smoothness scale rj decreases with Re. As 
the energy input increases the viscous energy dissipation 
rate, i/(Vt;)^, remains of the order of i//t^. Therefore far 
above the transition the principal part of the energy is 
dissipated by the polymer relaxation. Then the viscous 



term in Eq. (1.13) can be neglected and we obtain 



(3.1) 



is the wave frequency and k is its wave vector. Therefore 
the group velocity of these waves is B which can be es- 
timated in accordance with Eq. (3.1) as \/£t. The wave 
velocity fluctuates, but the fluctuations occur at times 
~ T and are slower than the wave oscillations at kr^ 3> 1, 
showing that the waves are well-defined. There exist two 
mechanisms of the elastic waves attenuation: polymer 
relaxation and viscous dissipation. The first mechanism 
leads to the scale- independent attenuation ^ t^^, which 
is smaller than the frequency, at fcr* ^ 1. The sec- 
ond mechanism produces the attenuation ^ i/k^, which 
is much smaller than the frequency for kr]^ <^ 1 where 
77* — iy{eT)~^/^. Thus the elastic waves attenuate weakly 
in the interval f;^^ < fc < 77^^. This interval can be called 
the elastic waves range. 

The dynamics at scales r ^ is also characterized by 
the stretching that takes place at a time scale r and is 
slower than the wave's oscillations. It is this dynamics 
that determines the velocity spectrum at fcr* 3> 1, since 
the wave oscillations do not influence the spectrum, like 
in the example presented at the end of Section As a 
result, we come to a power spectrum, which is examined 
in the next subsection. 



A. Power-law Spectrum in the elastic waves range 



where e is the energy injection rate per unit mass, esti- 
mated as V^/L. The relation (3.1) means that a typical 
value of B grows as the energy input increases, and, con- 
sequently, the elastic stress tensor does. 

The above quantities can be estimated, using the 
Kolmogorov-like reasoning. Then we obtain from Eq. 
(|t]) that B ~ Ver. Next, as follows from Eq. ( [1.12D , 
at the scale we have v ^ B. Equating then the 
characteristic velocity gradient v/r^ to l/r, one obtains 
~ (er^y^^. Note that this agrees with the direct 
Kolmogorov theory estimate of based on (5u(r*) ~ 
{er^y^^ ^ r^jr. Near the coil-stretch transition the 
viscous and ela stic d issipation terms in the energy bal- 
ance equation (1.13) are of the same order. Estimat- 
ing e by the viscous dissipation term vjr^ one finds 
RBc ^ [L/^ / {ut)]^/^ for the value of the Reynolds num- 
ber at the transition. The same answer can be found by 
equating r* and the Kolmogorov-41 estimate {u^/eY^'^ 
for the viscous scale 77. 

Below we investigate the case Re ^ Rcc that elucidates 
most clearly the role of the polymer elasticity. Since the 
condition implies that the viscous term is negligible at 
the scale r*, a new interval of scales, where viscosity is 
negligible but elasticity is not, has to exist below r^. The 
analogy with the magnetic hydrodynamics, noted above, 
helps us to understand dynamics of fluctuations in this 
interval. These small-scale fluctuations, which occur on 
the background of stresses excited at r ^ r.^, are elas- 
tic waves similar to the Alfven waves propagating in the 
presence of a large-scale magnetic field in plasma [^,^ . 
The dispersion relation for the waves is w = Bfc, where lo 



As we explained, statistical stationarity implies, that 
the velocity gradients fluctuate near l/r, the value char- 
acteristic of the scales r ^ r^,. Therefore velocity gradi- 
ents have to decrease with diminishing scale r at r < , 
which can be formulated as Vw' ^ Vw. Here v' is the 
small-scale component of the velocity containing only 
harmonics with wavevectors satisfying /cr* ^ 1 . The ex- 
istence of elastic waves at these scales leads to equipar- 
tition of kinetic and elastic energies (see and the 
proof below) so that VB' ^ VB holds too. As a result 
the influence of v' and B' on motions at scales ~ r* is 
negligible, that is v' and B' can be treated as passively 
advected and stretched by v and B. Equations for v' 



and B' can be found by linearizing Eqs. ( 1.1C , 1.12 ) 



dtB' + {v^)B' + {v'^)B = {B\/)v' + {B'\/)v - B'/t, 
dtv' + {vVy + {v'V)v (3.2) 
= -Vp' + {BV)B' + {B'V)B + vV^v', 

where p' — P' j q. 

The inequalities Vu' < Vt;, VB' < VB imply that 
at r < the differences 8v and bB, taken at points, 
separated less, than r*, scale linearly with the separation 
according to 



8vi 



5Bi = lijSvj 



(3.3) 



which is a generalization of Eq. (2.2). Both matrices 
<yik = VfeWi and 7^^ = VkBi have typical values l/r 
and correlation times of the order of r. To investigate 
the statistics of the small-scale components we may use 



7 



the same scheme, as was developed in Section || for the 
passive scalar, expanding the eq uat ions ( [3.2| ) near a La- 
grangian trajectory, like in Eq. ( ^.3| ). As we explained in 
Section ||, the velocity should be expanded to the first 
order (because of the Galilean invariance). The need to 
expand B to the first order too follows from the fact 
that zcroth order term produces Alfven waves that do 
not affect energy balance of the waves of the same type 
as explained at the end of Section ^ Passing to Fourier 
components (of the functions of the argument r ~ x), we 
obtain the equations 

= -ikp'i, + i[B ■ k)B'k + iB'k - yk'v'k, (3.4) 



i{B-kyk + <jB'k~—, 

T 



(3.5) 



analogo us to Eq. (2.4). The quantities cr, 7 and B in Eqs. 
( I3.4y3.5l ) are measured in the Lagrangian frame (they are 
functions of time and the Lagrangian marker). Corre- 
lation functions of the fields v' and B' are defined as 
averages over volume (or, what is the same, over differ- 
ent Lagrangian trajectories), that is over a statistics of cr, 
7 and B. There is a new ingredient in comparison with 
the consideration of Section ||, which is the stretching 
terms like av'k- However, these terms preserve zero scal- 
ing dimension of the time-evolution operator and, con- 
sequently, they are not expected to destroy the power 
character of the velocity spectrum E{k). 

From now on we neglect the viscous term in ( ^3.4|) which 
is justifiable for not too large wavevectors (a criterion is 
determined below). Using the incomprcssibility condi- 
tion, we express the pressure pj, — 2i[kav' — kjB']/k'^. 
The description is significantly simplified in terms of the 
Elsasser variables g± = v'k ± B'k, which satisfy 



dtg± - 
k fk 



k<y± • ^) 9± = ±*(S • f')9± - ^ (3.6) 



k \ k 



k fk 



where tT± = tr =F 7. Of course, the equations (3.6) are 
compatible with the conditions k ■ g± — following from 
the solenoidality of v and B . The right-hand side of ( |3.6| ) 
contains the Alfven term iB ■ k, describing the wave os- 
cillations. As we explained at the end of Section ||, it is 
convenient to eliminate the oscillations, introducing the 
corresponding phase and amplitude g± — a± exp[i(^±] 
with ip± satisfying 

dt^±- (^ka±--^^^±^±{B-k), (3.7) 

ip±=± f dt'B{t')W±^''^{t')W^{t)k, (3.8) 
Jo 



where dtW± = cr±W±, 
the amplitudes a± are 



W^±(0) = 1. The equations for 



df.a 



d \ k (k(T±a±) a± 



2t 



k{kazpa^) 



fc2 



exp(Tj<? 



(3.9) 



(/) = / dt'fc [W+{t)W^\t') + W-{t)WZ\t')] B{t') . 
Jo 

The above equations are, again, compatible with the con- 
straints k ■ a± = 0. 

We observe that a characteristic time of the variations 
of the amplitudes is r while the phase in the last term 
varies by 27r during the characteristic time {kB)~^ <^ r 
(the last inequality coincides with the previously derived 
condition for the existence of waves). Inde ed, the expo- 
nent appearing in the last line of Eq. (3.9) at t ^ r can 



be estimated as i3(t)Tfce*/'^, where fee /'^ is the current 
value of the wave vector, increasing due to the stretch- 
ing process. Averaging ( |3.9D over times much larger than 
{kB)~^ but much smaller than r (the procedure is noth- 
ing but the Bogolubov-Krylov averaging method), we 
find the following amplitude equation 



dta± 



d \ k f k \ a± ^ 



We observe that the equations for and a_ decouple. 
This is in accordance with the qualitative considerations 
of Kraichnan |2^ who argued that the interaction of the 
waves, described by the amplitudes and a_, is weak 
because their propagation directions are reverse. 

The equation (3.10) has the same structure as the 
equation for the linearly decaying scalar, considered in 
Section The difference is in its vectorial nature and in 
the presence of the term directed along k that comes from 
the solenoidality cond ition ka± = 0. A formal solution 
of the equation (3.10) can be written as 



a± {t, fe)=e-*/2-M± [t, q±]a± [0, q±] 
9tM±(t,fc) = (/^2/^) (/±a±Af±)^ 



(3.11) 
(3.12) 



where = kW±{t), f± — kW^^{t), and the ini- 
tial condition for the matrices M± is Mik{t = 0, fc) = 
Sik — kikk/k"^. The term kik^/k'^ in the initial condition 
for M± vanishes after contraction with solenoidal field 
a(t, fc), that leads to the correct initial condition for a±. 
Note that f±M±{t,k) = 0. Indeed, fcM±fO,fc) = and 
dt[f±M±{t,k)] = 0, as follows from Eq. j3.12| ) and the 
equations dtf± — ~f±<J±- 

Remind that B is defined up to sign. Therefore all 
the statistical properties of the solution have to be in- 
variant under the transformation B ^ B. This trans- 
formation interchanges g+ and g_. Therefore statistical 
properties of g+ and g_ are identical. Particularly, the 



8 



spectra of g+ and g- coincide, that is {g±i{k)g±j{k')) = 
{2TT)^6{k+k'){6ij-kikj/k'^)E{k). Thus, without any loss 
of generality, one can consider g^ solely. At calculating 
the correlation function of entering the definition of 
E, we can exploit an independence of W+{t) and 
of a+(0, k) that holds at t much larger than the correla- 
tion time T of cr+ (this is again in complete analogy with 
the consideration of Section Then the stationarity 
condition for E{k) reads 

2E{k)=cM-t/r) {Z{t,q+)E{q+)) , (3.13) 

kMl{t,k)A'I+{t,k)k 



Z{t, fc) =tr Mj^ {t, k)M+{t, fe)- 



fc2 



where we used M{t -k) = M{t,k). Note that Z = 2. 
Indeed, using Eq. (Ill) and f±M±{t,k) = one eas- 
ily shows that the time derivative of M'^{t, k)M{t, k) 
vanishes so that M'^ (t, k)M(t, k) — Sa — kikj/k'^. and 



Z(t,k) = 2. Thus the equation ( 3.13 ) simplifies to 

E{k) = eM-t/r) {Eiq+)) , (3.14) 

almost identical to Eq. ( |2.8|) established in Section 
||. Similarly, one can formulate an equation for the 
cross-correlation spectrum function E'{k) defined by 
{{g+Uk){g-)j{k')) = 5{k + k'){5^j-hkj/e)E'{k). The 
fast oscillating phase does not cancel in this equation 
leading to the inequality E'ik) ^ E{k). It leads to the 
conclusion that the spectra of v and B coincide and are 
equal to E{k)/2 each. This proves the equipartion of the 
elastic and the kinetic energies claimed above. 

In analogy with the consideratio n of Section || one 
can establish that the solution of (3.14) is power- like 
E{k) cx fc~" where a is determined implicitly by 



1 = exp(— t/r) 



|fcW+(t)/fc|« 



(3.15) 



Again, a > 3 (see Appendix |A 3| ). In fact a stronger in- 
equality follows from the stationarity condition. Namely, 
the spectrum has to decay faster than fc~^ {k^^ in 
the spherical normalization). Otherwise ((Vt>)^) = 
J E{k)k'^dk is determined by scales smaller than r*, vi- 
olating the condition that the gradients have to be satu- 
rated at the value 1 /t reached at . The condition a > 5 
coincides with the applicability condition of the above 
consideration that uses v{x + r) — v{x) « or for r <C r*. 
Indeed, then / E(k)k^dk is determined by fcr* < 1 and 
({v{x + r) — v{r) — crr)^) ^ {{ar)^) for r ^ r*. 

It is natural to ask, whether a is a universal num- 
ber, independent of Re. The above analysis shows that 
a is determined by the statistics of a and 7. Within 
the framework of the Kolmogorov theory the statistics 
is independent of the inertial interval length and can be 
characterized by a single parameter Ai — 1/t. We con- 
clude that the dimensionless quantity a is a universal 
number in the Kolmogorov theory. In fact, due to in- 
termittency the statistics of velocity gradients depends 



on the length of the inertial interval and, consequently, 
on the Reynolds number. The current understanding of 
intermittency does not allow us to estimate a for a given 
Re, yet some qualitative assertions can be formulated. 
Intermittency enhances the probability of large gradients 
which leads to faster transfer of energy to large wavevec- 
tors. Therefore the velocity spectrum becomes flatter 
as the intermittency increases. As we established, the 
length of the inertial interval decreases with the growth 
of the Reynolds number for polymer solutions (since the 
lower boundary of the inertial interval increases) . Con- 
sequently, the intermittency decreases as Re grows. We 
conclude that a should be a monotonically increasing 
function of Re. 

Now we establish the region of scales where the power 
spectrum exists. Its lower boundary is related to the vis- 
cous dissipation, which grows with increasing k, destroy- 
ing the power spectrum at large wavenumbers. Co mpar - 
ing the viscous term vk^v', the last one in Eq. (^.4|), 
with, say, the stretching term av' , we find, that the vis- 
cous term wins at the scale ~ y/in-- As a result the 
power-law terminates at fc '--^ (i^t)~^/^. For larger wave 
vectors the velocity spectrum diminishes faster than a 
power of k, that is the power spectrum occurs in the in- 
terval r;^^ < fc < (i^r)"i/2. Note that at Re > Rcc we 
have ^/lyr 3> ?7* so that the power spectrum occurs in the 
interval, where the elastic waves are well-defined. 



IV. ELASTIC TURBULENCE 

We pass to the case of low Reynolds numbers. Then 
a random (chaotic) flow can be excited due to elastic in- 
stabilities, if the Weissenberg number Wi is larger than 
unity. This is the situation of the recently discovered 
"elastic turbulence" 

We investigate the case Re ^ 1 where the substan- 
tial deri vative in Eg. ( 1.12| ) can be neglected. Then the 
system (L10,^tJ) becomes 

p-'^VP = {B\7)B + i^\7^v, \7-v = 0, (4.1) 
dtB + {vV)B = {BV)v - B/t, V • B = 0. (4.2) 



The inequality Re <C 1 implies, that the kinetic energy 
of the solution can be neglected in comparison with the 
elastic one. The dissipation of the elastic energy is, how- 
ever, due to both energy dissipation mechanisms (solvent 
viscosity and polymer relaxation): 

A j dr^ = -i j dr^-,. I dr{V,vu)\ (4.3) 



as fo llows from Eqs. (4.1,1^) and is in accordance with 
Eq. ( pri3| ). 

The system of equations (4T 4^ ) has to be comple- 
mented by the boundary conditions for the velocity, 
which in the absence of polymers would lead to Wi > 1 
so that the equilibrium state of polymers is unstable. 
The instability eventually leads to a chaotic, statistically 
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steady st ate m aintained by the non-linear dynamics of B, 
see Eq. (4^). Stationarity of the statistics, again, im- 
plies Ai = 1/t and stat iona rity of the pi — t/r statistics, 
as it stems from Eq. (4.2). It follows that the velocity 
gradients arc of the order of 1/t in the bulk. Therefore 
a boundary layer has to be formed, where the velocity 
gradient, exceeding 1/t, at the boundary, drops to the 
value 1/t in the bulk. In the boundary layer the flow 
is mainly shear and the polymers are weakly stretched. 
The existence of the boundary layer is observed also ex- 
perimentally [ p7[ . 

The above picture means, that instabilities lead to ve- 
locity fluctuations with scales determined by the size of 
the system and a characteristic gradient 1/t. The fol- 
lowing estimates for the values of these large-scale fluc- 
tuations hold 



L 



L 



T 



— ' 152 Re«l, (4.4) 

i}^ VT 



where L is the linear system size. The estimates (4.4) 
imply that both dissipative terms in Eq. (^^) are of the 
same order. The correlation time of B and v is deter- 
mined by the typical value of the stretching and is of 
the order of r. The large-scale velocity fluctuations pro- 
duce smaller scale fluctuations of B that in turn induce 
small-scale fluctuations of velocity. The velocity gradi- 
ents become smaller when the scale decreases, since the 
large-scale velocity gradient is of the order of 1/r and 
the total velocity gradient is fixed at this value by the 
stationarity condition. This is in complete analogy with 
the consideration of the previous section. 

Now we are going to consider statistical properties of 
the small-scale fields v' and B' . The fields evolve pas- 
sively in the large-scale fields v and B. However, there 
is a major qualitative difference from the high-Reynolds 
case which is in the role of the large-scale component of 
B. In the high- Reynolds number case (see Section III) 
the terms with the large-scale field B in the equations for 
the small-scale fields conserve energy and, c onseq uently, 
do not enter the equation for the spectrum ( ^.13 ). That 
is why the correct description of the dynamics required 
account of the small gradient V^B on the background 
of B. In the elastic turbulence case the terms with B 
are dissipative and, consequently, one can forget about 
the gradient terms, as subleading ones. Thus the equa- 
tions for the s mall-scal e components of the fields, follow- 
ing from Eqs. (4.1,4.2), are 



Vp' = (BV)B' + i^V^v', (4.5) 
dtB' + {vV)B' = {B'W)v + {BV)v' - B' /t . (4.6) 

Let us stress that the above equations assume only B 
B' (not Vi? ^ VB' , which is in fact wrong for the elastic 
turbulence) . 

To analyze the above equations for the small-scale com- 
ponents we use the same scheme as in the previous sec- 
tions. As we explained above there is no need to account 



for the spatial variation of B so t hat (for the Fourier 
components) the equations (4^,4^) take the form 



, i(k-B). 

v' = ^ -B, p' = 0, 



dtB'~ 



B' {k-Bf ^,^ 



(4.7) 
(4.8) 



The constraints k ■ B' — k ■ v' ^ that follow from 
V ■ t) = = V • B are consistent with the equations 
(4.7 4.J). The last term in the equation for B' is due to 
the viscous dissipation, the term is of the same order as 
the elastic dissipation term —B'/t as follows from ([1.4|). 
We observe that, again, the time-evolution operator for 
B' has zero scaling dimension and (in accordance with 
the discussion at the end of Section |l^ one expects that 
the spectrum F(k) of B obeys a power-law. To demon- 
strate this pow er behavior we use a formal solution of the 
equation ( [l.q ) 



B'{t,k) ^ WB'{Q,kW)eyiY> 



t 



/ dt'C(i') 



(4.9) 



vS,{t') = [B{t')n{t', t, k)f, dt'Ti = -ncT + n{nan), 

where W — W{t,0), and n{t',t,k) is determined by the 
above equation with the final condition nk(t, t, k) = k/k. 

Let us analyze the stress spectrum function F{k): 
{B'.{k)B'j{k')) = S{k + k')F{k){Sij - kikj/k'^). Assuming 
that t is much larger than the correlation time t of B 
and cr, we may average independently over the velocity 
statistics at negative and positive times, as we did in the 
previous sections. Then we find 



2F{k) = (Zexp 
Z = exp(-2t/r) 



-2 fdt'm 
Jo 



kWW^k 



F{kW)) , (4.10) 

Tl2i 



(4.11) 



The coefficient of F in the rig ht-hand side of Eq. ( |4.10| ) 
is ind epen dent of the absolute value k oi k. Therefore 
Eq. ( 4.10 ) admits a power solution F{k) cx fc~^, where 
f3 has to be determined from the equation 



2 = 



Z exp 



^2/Jdi'^(t') 



\kW{t)/kf 



(4.12) 



Note that at t ^ t the average in Eq. ( 4.12 ) is deter- 
mined by the events with fc directed along the eigen vec- 
tor of the matrix WW^ , corresponding to the smallest 
eigen value (in this case the denominator in the expres- 
sion achieves a minimu m). T hen the second term in the 
right-hand side of Eq. ( [4.1l| ) can be neglected, and we 
obtain Z « exp(2/9i — 2i/T) (therefore, due to the station- 
arity of the pi -^t/T probability distribution, a statistics 
of Z is time-independent). The positive noise ^(t) can 
be considered as stationary with the correlation time r. 
Indeed, it follows from Eq. (4.9) that n{t') forgets its 
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final direction (and thus also becomes fc-independent) at 
t — t' ^ r. Thus the situation is similar to the one 
analyzed in Sections P. p^II A . Again, we can prove the 



inequality /3 > 3, see Appendix A 3 



We are now in a position to establish the spectrum of 
the velocity E{k). Indeed, it follows from Eq. (4/7) that 
E{k) (X k~^~'^ and the spherically normalized spectra 
obey 

Esphik) - v^L{kL)-^, Fsphik) - B^L{kLf-^ . (4.13) 

The estimates for the values and are written in 
Eq. (4.4). We observe that our scheme that assumes 
Vu' ^ Vw and B' B is self-consistent due to the 
inequality /3 > 3 proved above. It agrees with the exper- 
iment, where (3 is 3.3 3.5 

Finally, let us discuss the question concerning the va- 
lidity region of the power spectrum. Probably, it is deter- 
mined by the finite diffusivity k of the polymer molecules 
which is described by adding the term KV^Hij to the 



right-hand side of Eq. (1.2). Comparing this term with, 
say, the relaxation term with r, one concludes that the 
power spectrum terminates at fc^if ~ Vkt. At smaller 
scales the velocity spectrum diminishes much faster due 
to diffusivity. 



V. CONCLUSION 

We have investigated properties of turbulence in dilute 
polymer solutions in the cases where polymer molecules 
are strongly stretched. We established power-law distri- 
butions of kinetic and elastic energies over scales in some 
regions, where these power-laws are not related to an 
energy or other conserved quantity cascade (in contrast 
to the usual turbulence). In fact, excitation of elastic 
degrees of freedom at any scale leads to energy dissi- 
pation since the elastic dissipation is scale-independent. 
However, precisely this scale-independence can lead to 
a scale-ivariance in the dissipative intervals, where the 
flow can be treated as smooth. Small-scale fluctuations 
are relatively weak and evolve passively in the smooth 
flow. As a result, the evolution of fine-scale fluctuations 
depends trivially on the scale and power-law spectra are 
formed. Let us now describe the cases where the above 
general ideas are applicable. 

The first case, we examined, is the high Reynolds num- 
ber flow above the coil-stretch transition, when elastic 
degrees of freedom are activated. Strong interaction be- 
tween the elasticity and the flow modifies the latter be- 
low the scale (at this scale the velocity gradients are 
of the order of 1 /t) , which is the new energy dissipation 
scale. This scale is of the order of the Kolmogorov scale 
at the transition and becomes larger as Re is increased. 
At r ^ the properties of turbulence are the same as in 
Newtonian fluids. The energy cascades downscales from 
the pumping scale and dissipates due to polymer relax- 



principal Lyapunov exponent Ai fixed at the value 1/r by 
the elastic back reaction. The fluctuations in the interval 
of scales z^(er)~^/^ ^ r ^ are elastic waves. This leads 
to the equipartition of the kinetic and of the elastic en- 
ergies, that is the velocity spectrum E{k) and the elastic 
spectrum F{k) coincide at these scales. The smoothness 
of the flow at r < leads to the conclusion that these 
spectra are power-like and in the spherical normalization 
decrease faster than k~^ at kr^, > 1. The power spectra 
terminates at the scale (i/r)^/^, where the viscous dissi- 
pation overcomes stretching. 

The second case, we examined, is the elastic turbu- 
lence regime 1^-^. It is a chaotic state which is realized 
at small Reynolds numbers Re. The velocity gradient im- 
posed on the system by the boundary conditions exceeds 
1 /r which activates polymer degrees of freedom leading 
to hydrodynamic instability and chaotization. Again, the 
power spectra F{k) and E{k) are power-like in this case. 
However there are no elastic waves that would lead to 
equipartition. The main energy is carried by the poly- 
mers: E{k) ~ Re {kL)~^F{k), where L is the size of the 
system. The velocity spectrum E{k) decays faster than 
k~^ in the spherical normalization, which corresponds to 
the experimental data 

The above-described mechanism of forming power-law 
spectrum for small-scale fluctuations in a chaotic flow 
seems rather general to be realized for other systems. We 
expect it to occur in certain regimes in magnetohydrody- 
namics, flows in liquid crystals and low-dimensional flows 
on a substrate. 
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APPENDIX A: LONG-TIME LAGRANGIAN 
STATISTICS 

Let us briefly review the long-time statistical proper- 
ties of the Lagrangian mapping matrix W, determined by 
Eqs. (1.5 1.6). We consider W{t^,t_) at t+ > t_ and as- 



sume that tj^ — t- is much larger than the Lagrangian cor- 
relation time To- of the velocity derivatives matrix ( |l.(j ), 
where one expects a universal statistics . If the veloc- 
ity statistics is homogeneous in time, the probability dis- 
tribution of W{t ^^t^ ) depends on the difference t+ — t_ 
only. Equation (1.5) implies that at t+ — t- ^ the 
matrix is a product of a large number of independent 
matrices. This is the main reason for the universality of 
the W statistics. 

It is convenient to decompose the matrix W as 



ation at r 



The flow is smooth at r < with the 



W{t+,t^) = MAN, 



(Al) 
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where A is a d iag onal matrix, and M and N are orthog- 
onal matrices |2^. We denote the diagonal elements of 
A as e''^ , e''^ , and e''^ , and assume that they are or- 
dered: pi > P2 > P3- As a consequence of the c onst raint 
det W = 1 we have pi -f P2 + Pa = 0. Equation (1.5) can 
be rewritten in terms of pi, and the matrices M and N. 
The equations for pi are 



dpi/dt+ = di, 



(A2) 



where a = M^aM and no summation over the repeat- 
ing index i is implied. The matrices M and N satisfy 
dtN = riiN and dtM = MrJa, where 



1 )ik 



2sinh(pj - pk) 



o2pk 



p2pi 



iov i ^ k and (rii)ifc = (1^2)ife = for i = fc. It is possi- 
ble to show that the eigenvalues of W repel each other, 
so that the inequalities e'^i ^ e''^ ^ e''-'* are satisfied at 
t+ — t_ ^ To- 30|. Then the matrix ili tends to zero 
exponentially fast, i.e. N is determined by times of the 
order of in the vicinity of The matrix Q2 becomes 
/9-independent at i+ — t_ ^ Tq- and the evolution of M 
is decoupled from that of pi. Then the value of M is 
determined by the time of the order of at t w i+, i.e. 
at — t- ^ it becomes t_ -independent. 
The solution of Eq. (^) is 



dt' au{t'). 



(A3) 



where the right-hand side of Eq. (A3) is an integral of a 
random process independent of pi. Equation ( |A3D shows 
that the variables pi fluctuate around their average val- 
ues Ai(t+ — t-). Here the constants are equal to (an). 
They are called the Lyapunov exponents of the flow. 
Generally, the spectrum of the Lyapunov exponents is 
non-degenerate: Ai > A2 > A3, which is a necessary con- 
dition for the formalism to be self-consistent. The incom- 
pressibility condition ensures the identity A1-I-A2-I-A3 = , 



which implies Ai > and A3 < 0. Using the relation ( |l.7| ) 
one can show that Ai is indeed the average logarithmic 
divergence rate of two nearby Lagrangian trajectories: 

(din I (5a; I /di) = Ai . 



Lagrangian Statistics in usual Turbulent or 
Random Flows 



Here we consider the Lagrangian statistics in the case 
of a usual random (turbulent or chaotic) flow, when the 
velocity has finite correlation time and no constraints are 
imposed on the flow. Then the quantity pi can be treated 
as a sum of a large number of independent random vari- 
ables, provided — ^ t„. It is known from the sta- 
tistical mechanics (see, e.g., |3l)) that the distribution of 
such quantities is given by the exponent of an extensive 
function. In our case the probability distribution func- 
tion (PDF) of p, is 



'P{t,Pi,P2,P3) oc iexp 
xS{pi + P2+ ps) , 



-tS 



Pi - Alt p3 - Xst 



t 



t 



(A4) 



where t = — and pi > p2 > Ps is implied pOji . The 
main exponential factor of the PDF has a self-similar 
form described by the function S, which can be called 
entropy function (see 32 3^). It is positive, convex 
and has a minimum at zero values of its arguments. The 
precise form of S is determined by details of the velocity 
statistics. The PDF has a sharp maximum at pi = Xit. 
In its vicinity the function S has a quadratic expansion, 
i.e. the distribution of p is Gaussian. However, if one is 
interested in the expectation values of exponential func- 
tions of Pi, they are determined by the wings of the PDF 
where the Gaussian approximation is invalid. This en- 
tails the use of the whole entropy function. 

To average the functions of pi only, one can introduce 
the reduced probability distribution function. 



V{t,pi) cx — exp 



-tSi 



Pi - Alt 
t 



(A5) 



which is an integral of V{t, pi, p2, ps) over p2 and ps. At 
small X the function Si{x) can be written as 



(A6) 



Here A — Jdt (((Th (t)(Tii (0))) (where double brackets 
designate irreducible correlation function) determines the 
dispersion of pi: ((pi — Ait)^) ~ tA. Expansion (A6) is 
sufficient to describe typical fluctuations of pi, whereas 
the whole function S is needed to describe rare events. 



Similarly, Ai -I- A2 = — A3 is the average logarithmic rate 
of the area growth. 

Note that at t+ — t_ ^ the statistics of M, A and 
N are independent. Indeed, the values of pi are accumu- 
lated during the whole evolution time <+ — t_ (see ([A^)) 
and are not sensitive both to the interval (t-,t- + Tq.) 
determining N and interval (t+ — r^, t+) determining AI. 
Both matrices M and N are distributed isotropically be- 
cause of the assumed isotropy of the velocity statistics 
(this is the isotropization of W referred in the main text). 



2. Special properties of the long-time Lagrangian 
statistics above the coil-stretch transition 

Above the coil-stretch transition the Lagrangian statis- 
tics acquires new qualitative features caused by the poly- 
mer back reaction. They can be inferred from the equa- 
tion ( |l.9| ) which leads to 

In \Bit, x)\ - In \Bito, r)\ « pi(t, to) - (A7) 
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The expression (A7) is correct, provided t — to^^r (since 
the polymer relaxation time r determines also the ve- 
locity gradients correlation time). The left-hand side of 
the equation (A7) has stationary statistics which leads 
to dramatic changes in the statistics of pi . Upon averag- 
ing one finds (pi) = t/r which means Ai — (cth) — 1/t 
as it was already explained in the main text. Stationar- 
ity of the dispersion of pi leads to the conclusion that 

di'((o-ii(0)(Tii(t'))> fx t-^ in the limit t oo. This 
is related to the anticorrclation property of an, men- 
tioned in the main text. That means that the dispersion 



A, defined by Eq. (A6), vanishes above the coil-stretch 
transition. In fact, this vanishing is not abrupt and oc- 
curs within 1/ ln(i?back/-Ro) vicinity of Re^. Finally one 
concludes that ((tTii(O)CTii(i))) oc t^"^ at large t which is 
again very different from a Newtonian fluid where expo- 
nential decay is observed, see Ref . S . 



3. Inequality for the exponents a, [3 

Here we establish the inequality for the exponent a, 
characterizing the passive scalar spectrum, see Section 
U, and the inequality for the exponent (3 appearing in 
the elastic turbulence problem, see Section IV. In both 
cases we investigate the solution A of the equation of the 
type 



IfeW^I'^exp 



dt'y{t') 



i A 



(A8) 



where t is much larger than the correlation time of a 
random positive noise y{t) and a. At these times the 
behavior of the moments is exponential and one can de- 
fine {eyi^[- f^dt'y{t')]\kW\^) ~ fc"^ exp[7(5)t] where ^{5) 
is a convex function due to Holder inequality. This func- 
tion is strictly smaller than another convex function 7 ((5) 
defined by {\kW\^) - fc* exp[7(5)t]. 

The function ^{5) has a universal behavior which we 
describe now. It is convenient to write 



{\kW\^)^ / dk'k'^ {5{k~k'W-\t))) , 



(A9) 



making it explicit that the wavevectors evolve accord- 
ing to k{t) = k{0)W-\t) ||. Introducing k'W~^{t) = 
k' exp [p{t)] n, where n is a unit vector, one finds 



Pit) 



dt'at') 



An 
dF 



= na + nC, Q = nan. 



One observes that C is independent of p{0) which leads 
to the conclusion that at t 3> ''"c(o'), where Tc{(j) is the 
correlation time of cr, the probability distribution of p is 
described by an entropy function S. 

It can be shown that (p(t)) = jAsji, where \\^\ is the 
lowest in the hierarchy of the Lyapunov exponents of the 
fiow. This fact is intuitively clear as A3 determines con- 
traction in the real space and thus stretching in fc-space. 



We note that p{t) is determined by the whole interval 
(0,t) while n{t) only by Aj~^ vicinity of t. As a re- 
sult, p and n are independent at Xit 3> 1. Since n is 
isotropically distributed over the unit sphere one finds 
{6[k^k'W-^{t)]) ^ (5T£-ln(A:/A:'))exp[-3p])/(47rA:'3). 
Substituting this into (A9) and performing the integral 
one finds 



dp 



-{5 + i)p~tS 



t 



exp[7((5)t] 



where N oc \ft is the normalization factor insignificant 
for the following considerations. It follows that besides 
the zero at the origin following from the definition, the 
function 7((5) vanishes at (5 = —3. This is a general con- 
sequence of the isotropy employed above, see also 
Besides, we observe that 7'(— 3) = A3 < 0. Combining 
these properties with the convexity of 7(5) we conclude 
that 7(5) is negative for —3 < (5 < and positive other- 
wise. 

To prove the inequality on A appearing in Eq. ( [A^ ) 
one notes that both 7(5) and 7(5) tend to 00 as \b\ 00. 
Then it follows from j{6) < ^(6) that there are two solu- 
tions of Eq. (|A^): one positive and one smaller than —3. 
Substituting y{t) ~ 2 /t and recognizing that the expo- 
nent a appearing in (p^) must be positive we conclude 
that a > 3. Analogously, substit uting y{t) = 2^(i) we 
conclude that the solution of Eq. (4.12) satisfies /3 > 3. 

Finally, let us give an example of calculating a in a 
limiting case. The exponent is determined by the equa- 
tion 



dp 



{a — 3)p — tS 







'2t 




= exp 






T 



(AlO) 



Note that a — 3 vanishes in the limit \X3\T — s- 00 sin ce in 
this limit the linear decay term is negligible in Eq. (2.1) 
and Batchelor k^^ spectrum must result. Therefore at 



large IA3IT the integral (AlO) is determined by the max- 
imum of the probability concentrated at p = \X3\t which 
leads to 



-3-2(|A3|r)'- 



|A3|t> 1. 



(All) 



For a general value of lAsjr the exponent a is determined 
by the concrete form of the entropy function. 
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